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Abstract. Magnetic field modifies the properties of waves in a complex way. 
Significant advances has been made recently in our understanding of the physics 
of waves in solar active regions with the help of analytical theories, numerical 
simulations, as well as hi-resolution observations. In this contribution we review 
the current ideas in the field, with the emphasis on theoretical models of waves 
in sunspots. 



1. Motivation and questions 

Photospheres and chromospheres of sunspots are the regions where different 
physical agents come into play with almost equal weight. The restoring forces 
in the wave equation, such as magnetic Lorentz force, gas pressure gradient and 
buoyancy, are of the same order of magnitude, making arbitrary the division 
into pure wave modes and complicating theoretical models. A recent summary 
on the observational properti es of sunspot oscil l ations together with some the- 
oretical issues is presented by iBogdan Sz Judge] (|2006l ) . The most recent obser- 
vational results suggest that wave phenomena at different lay ers of the sunspots 
umbrae and penumbrae are related with each other ( see e.g. iMarco et al] | l996 |; 
Maltbv et al.lll999l.l200H:lBrvnildsen et al.ll200dl2002l:IChristopoulou et al-lboocl. 
20011: iRouppe van der Voort et al.ll2003l : iTziotziou et al.1 12006 : iBloomfield et al.l 



2007al ) . The interpretation of the observational material defines several groups 
of questions concerning the sunspot wave physics that can be summarized in the 
following way: 

• What drives the waves observed in sunspots? Are they externally driven 
by the quiet Sun p-modes? Are there sources of oscillations inside the 
umbra due to weak convection? 

• How can the oscillations observed at different sunspot layers be interpreted 
in terms of MHD waves? What are the relationships between photospheric 
and chromospheric oscillations? What causes the complex spatial pattern 
of oscillations such as chromospheric umbral flashes, penumbral waves, 
spatial coherency of waves over the umbra, etc.? 

• Why is the wave power suppressed in the umbral photosphere compared 
to the quiet Sun? Why is it enhanced in the chromosphere? 

• What mechanisms produce the change of the dominating frequency of 
waves in the umbra from 3 mHz in the photosphere to 5-6 mHz in the 
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chromosphere? What are the reasons for the spatial power distribution at 
different frequency intervals observed in the umbra and penumbra? 

• Are there observational evidences of the mode transformation in sunspots? 

• What is the source of helioseismological velocity signal detected in active 
regions? Are they due to fast, slow MHD waves? What are the con- 
sequences of the strong magnetic field of sunspots onto helioseismology 
measurements and inversions of sub-photospheric structure of sunspots? 

Some of these questions, like the interpretation of the wave modes or change of 
the frequency with height, received larger attention in the literature in the past 
and can be considered answered to a smaller or larger extent. Some others, like 
those concerning local helioseismology in active regions, only received attention 
recently and no unique answer has been found yet. In the rest of paper we will 
describe the existing models and the theoretical progress made in understanding 
sunspot oscillations. 



2. Excitation of oscillations in sunspots 



It is generally accepted that the quiet-Sun oscillations are s tochastically excited 
by so l ar convection at the top part of the convection zone (jGoldreich Kumar 
1990l : iNordlund Sz Steinl l200ll ). The frequency spectrum and the temporal be- 
haviour of waves in the umbra is ver y similar to the quiet Sun, except for reduced 
power (e.g. iPenn &; LaBontd fl993l ) . It is systematically found in observations 
that more wave power travels toward the umbra than leaves the umbra. One of 
the natural explanation of this effect is to assume that the sunspot waves are 
driven by external p-modes, modified inside the magnetized sunspot atmosphere. 

A number of theoretical investigations has been performed assuming the 
incom ing p-mode waves incident on a sunspot-l i ke magnetic field concentr ation 
(e.g. ICallv et al.l 1 19941 : ICallv & Bogdanl Il997l : iRosenthal & Julienl l2000l ) and 
studying the wave transformation, "absorption" and scattering by sunspots. 
These studies have shown that both amplitudes and phases of the incoming 
waves are modified after passing thought the magnetic field regions. In particu- 
lar, the surface amplitudes of the p-modes were found to be reduced due to the 
mode transformation (Sect. 5). 

Despite the presence of strong magnetic field in sunspot's umbra and penum- 
bra, the m odels of magnetoconvection indicate that these environments can be 
connective ( Ledll993l; Weiss et al1ll990l : ICattaneo et al.1120031 ; ISchiissler fc Vogler 
2006j; ISpruit fc Scharmeri l2006f T The narrow nearly field-free upflowing umbral 
dots adjacent by downflows are the result of convection in a str ong magnetic 
field. Recent observations seem to confirm this idea (see e.g. iBharti et al 



2007; IWatanabe eFall I2008T). Thus, so me w aves can be generated inside the 
sunspot umbra as well. Jacoutot et al. ( 20081 ) performed numerical simulations 
of magnetoconvection and studied the spectra of generated waves as a function 
of magnetic flux in the model. Apart from the power suppression in regions with 
enhanced magnetic field, these simulations suggest an increase of high-frequency 
power (above 5 mHz) for intermediate magnetic field strengths (of the order of 
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300-600 G) caused by changes of the spatial-temporal spectrum of turbulent 
convection in a magnetic field. 

It can be concluded that both external driving and in-situ excitation by 
convection are able explain the observed properties (e.g. power reduction) of 
waves in sunspots. A qualitative estimation of the effect of the reduced wave 
ex citation in sunspots and th e direc t comparison to observations was reported 
by iParchevskv Sz Kosovichevl ( 20071 ). The authors performed hydrodynamical 
simulations of waves generated by random sources, with the strength of the 
sources reduced to zero in the sunspot umbra. They obtained that, even though 
no waves were excited in the umbra, the velocity power measured there was about 
twice larger than actually observed. Thus, the waves detected in sunspots must 
be a mixture of external p-modes (and other wave types produced by the mode 
conversion) and MHD waves generated directly inside the sunspot by the weak 
convection. The relative contribution of the different effects is to be determined. 



3. Interpretation of oscillations in terms of MHD waves 



Historically, the sunspot oscillations were divided into 5-min photosph eric oscil- 

lations, 3-min chromospheric oscillations and running penumbral waves ( Bogdan fe Judg"el 
2006). Recent observations show that these oscillations can be a manifestation 



of the same physical phenomena (IMaltbv et al.l 11999. 120011; iBrv nildse n et al 



20001. 120021: IChristopoulou et aUl200CL l200ll: IRouppe van der Voort et alJ l2003b 



Centeno et alJ 12006b iTziotziou et al.ll2006l . 120071 ) Chromospheric observations 



show umbral flashes followed by a smooth continuous quasi-circular expan- 
sion of the visible perturbation pattern from the umbra through the penum- 
bra toward the edge of sunspot. At the same time, the frequency of chromo- 
spheric oscillations dec reases from 5 to 3 mHz from the umbra to the penumbra 
(|Tziotziou et al.1 120071 ). The waves in the umbra propagate from the photo- 
sphere to the chromosphere along the magnetic field lines forming shocks of 
about 10 km/sec peak-to-peak amplitude at heights of formation of He I 10380 
A line (lLitedll986l: i Looez Ariste et~aD l200ll: IRouppe van der Voort et alJl2003l : 
ICenteno et alJl200d ). 

The characteristic speeds of waves change many orders of magnitude along 
the photosphere and chromosphere of sunspots. At some height waves propagate 
through the region where the sound and the Alfven speed are equal (eg = va). 
In this region mode transformation and coupling of different wave phenomena 
occurs. In addition, the magnetic field and thermodynamic variables show im- 
portant gradients both in horizontal and vertical directions. All these ingredients 
make realistic modeling of sunspot waves a rather difficult task only accessible 
via numerical simulations. Theoretical bases of the MHD wave propagation in 
a stratified atmosphe r e wit h a constant magnet i c field were developed by, e.g. , 
Ferraro Plumptonl ([1951): lOsterbrock! (|l96lh : IZhuezhda fc Dzhalilovi (|l98l 
19841 ): lEvans Roberts! (jl99r£ ~ One of the first numerical simulations of waves 
in non-trivial magneto-atmosphe res with applicat i ons to solar photosphere an d 
chromosphere was performed bv iRosenthal et al.l (120021 ): iBoedan et all (|2003T ). 
Similar calculations for conditions appropriate for su nspots rather than flux 
tubes were carried out bv lKhomenko &: Collado"sl ( 20061 ). The waves were gener- 
ated by a photospheric pulse located inside the magnetic field region and propa- 
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gated through the layer where c $ = va on t heir way from t he ph otosphere to the 
chromosphere. iRosenthal et al.l ( 20021 ) and iBogdan et all ( 20031 ) pointed out an 
important role of the magnetic canopy defined as the region where the plasma 
/?(/? = Sirp/B 2 = ( / y/2)c 2 s /v\) is close to unity. According to their results, the 
properties of waves observed in the atmosphere depend sensitively on the loca- 
tion and orientation of the magnetic canopy. The mode transformation occurs 
in the canopy region and several wave types can be present simultaneously there 
traveling in different directions. However, the fast and the slow wave modes 
decouple effectively away from the transformation layer. At hei ghts above (5 = 1 



the r e fraction and reflection of the fast (magnetic) mode occur ([Rosenthal et al 



|2002| ; iKhomenko Colladod 20061 ). The slow mode (either produced by the 
driver or after the mode transformation) is always guided along the magnetic 
field lines, except at t he region /? ~ 1 where it deviates from this direction by 
maximum 27 degrees ( Osterbrocklll96l1 ). These and other theoretical investiga- 
tions demonstrated that the explanation of the observed oscillations in magnetic 
structures in terms of MHD waves requires the knowledge of the location of the 
j3 = 1 level, as well as the inclination of the magnetic field lines. 

Mathew et al. ( 20041 ) used spectropolarimetric observations in Fei 1.56 /im 



lines to calculate maps of plasma (3 for the observed sunspot. The plasma beta at 
heights of the continuum formation of radiation at 1.56 /xm was found to be less 
than 1 in most of the umbra and inner penumbra reaching values above 1.5 in 
the outer penumbra. Most photospheric lines used typically in observations (like 
Nil 6768 A used by SOHO/MDI) are formed even higher in the atmosphere. It 
indicates that oscillations observed in sunspot umbra should be due to the low (5 
MHD waves. The situation is more complex in the penumbra where the plasma 
(3 should be around unity at the typical heights of formation of photospheric 
spectral lines. 

From observations, there is no clear conclusions on what wave types are 
detected in sunspots. The analysis of magnetic field oscillations in a sunspot 
obtained from spectrop olarimetric observations of Fei 1.56 /jm lines allowed 
Khomenko et al.l (|2003l ) to conclude that the observed waves are a mixture of 



fast and slow low-/? waves and the contribution of the slow (acoustic-like) wave 
is the dominant in the sunspot umbra. The slow (acoustic-like) wave propaga- 
tion along the magnetic field lines from the photosphere to the chromosphere 
in the sunspot umb ra was reported fr om spectropolarimetric o bservations by 
Centeno et al.l ( 20061 ). In the penumbra. iBloomfield et all ( 2007bl ) obtained that 



the propagation speeds of waves in the deep photosphere most closely resemble 
the fast-mode speeds of the modified p- modes, i.e. high-/? acoustic-like waves 
propagating towar d the sunspot at an angl e ~ 50 degrees to the vertical. In a 
subsequent paper, IBloomfield et al.1 (|2007al ) used a time series of photospheric 
Si I and chromospheric Hei spectra at 1.08 /im and concluded that the running 
penumbral waves are a visible pattern of low-/? slow mode waves propagating and 
expanding their wavefront along the inclined magnetic field lines in the penum- 
br a. Similar interpre t ation of the running penumbral waves was also suggested 
by IBogdan Judge! ( 20061 ). Anticipating the conclusions from Sect. 6 of this 
paper, theoretical modeling of local helioseismology signals in sunspot regions 
seem to suggest that they are rather due to high-/? fast mode waves (modified 
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p-modes) ( Khomenko et al. 2008). Given these diverse conclusions, more work 
is required to unify the picture of MHD wave modes observed in sunspots. 



4. Frequency distribution 

The power of velocity and intensity oscillations is distributed in a complex way 
over the active regions dHindman & Brown|[l998l : I Jain Habe3 2002: T ziotziou et al 



20071 ; iNagashima et al.ll2007l ; lMathewll2008l ). In the photosphere, the power of os- 



cillations in generally suppressed at all frequencies, both in the umbra a nd in the 
penu mbra, except for the bright ring at the umbra-penumbra boundary (jMathew 
20081 ) . Acoustic halo s are found at the edges of active regions depending mag- 



netic field strength ( Braun et al. 19921 ; iBrown et al.l 119921 : iHindman Brownl 
1998) being most prominent at high frequencies above 5 mHz. In the chro- 
mosphere, the power is suppressed at low frequencies (below 2.5 mHz) in the 
umbra but then is en hanced in the high frequencies (above 3.5 mHz) (e.g. 
INagashima et aill2007l ). The most prominent are 5-6 mHz oscillations observed 
in the umbral chromosphere. The run ning penumbral waves in the chromosphere 
have lower frequency around 3 mHz ( Tziotziou et al.ll2007l ). 

Assuming the waves observed in the umbral photosphere and chromosphere 
are slow low-/? acoustic-like waves propagating along nearly vertical magnetic 
field lines, the mechanism of the frequency variation with height of these waves 
shoul d be similar to t h at of the ordinary acoustic waves in the non-magnetic 
Sun. Fleck Schmitz ( 199ll ) have demonstrated that the change of frequency 
with height from 3 to 5-6 mHz is due to a resonant excitation at the atmospheric 
cut-off fr equency, occurring ev en for linear waves in an isothermal atmosphere 
(see also iKalkofen et al. 19941 ). In the solar case, the low temperatures in the 
upper photosphere give rise to the cut-off frequency around 5 mHz. The re- 
sponse of the solar atmosphere to non-linear adiabatic shock wave propagation 
also leads to an appearance of the 5 mHz frequency peak in the power spectra, 
under t he condition the underl ying photosphere has a low frequency 3 mHz com - 
ponent (Fleck Schmitzlll993r). At the same time, IZhugzhda fc Loca"nsl ( 198ll ); 
Gurman fc Leibacheri (|l98^ : IZhugzhdal (|2007t ) argue that the observed spectrum 



of sunspot umbral oscillations in the chromosphere is rather due to a particular 
temperature gradients of the atmosphere acting as an interference filter for linear 
three-minute period aco ustic waves. Yet another explanation was suggested by 
Carlsson Steinl ( 19971 ) based on the fact that at lower frequencies the energy 



falls off exponentially with height, stronger that for the high -frequency waves 
not affected by the cut-off. Similar effect was discussed also bv lBogdan &: Judge 
(|2006ft . 

More puzzling is the horizontal distribution of wave frequencies over ac- 
tive regions observed at different heights. The power suppression in the um- 
bral photosphere is thought to be due to t he mode convers i on extracting the 
energy of the incident p -mode waves (e.g. ICallv fc Bogdanl Il997l : ICallv et al " 



2003; ICrouch et all 120051 . see the discussion in the next Section). The power 



distribution in the chromosphere can be explained assuming again the low-/3 
slow mode propagation along the inclined magnetic field lines (jBogdan &: Judge 
20061 ; iBloomfield et al.l l2007al ) . In the low-/? plasma, the effective cut-off fre- 



quency along each individual magnetic field line decreases by a cosine of its 
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inclination angle 9. Taking the typical values of the inclination in the umbra 
and penumbra, this mechanism is able to produce the required change from 
the 5 mHz frequency characteristic for umbral flashes to 2-3 mHz characteristic 
fo r running penumbral way es, as confirmed by spectropolarimetric observations 
of Bloomfield et all ( 2007al ). As for the high-frequency acoustic halo surround- 
ing active regions, n o plausible expla nation exists as of today. In his recently 
reported simulations iHanasogd (2008) calculated the rms wave power after ran- 
domly distributing acoustic sources over the region containing a sunspot-like flux 
tube. Both power reduction in the "umbra" and halo effect in the surrounding 
appeared naturally in these simulations, confirming that they are definitely an 
MHD effect. Certainly, the mode conversion and wave refraction in the inclined 
magnetic field at the edges of active regions play a role in the appearance of 
acoustic halo, requiring a further study. 



Mode transformation 



Mode transformation causes important effects on the observed wave propaga- 
tion in active regions. The bases of the MHD wave transformation theory i n 
a stratified solar atmosphere were developed bv IZhuezhda & Dzhalilovl (|1982T ). 
Approaching the layer where the acoustic and the Alfven speeds are equal, the 
phase speeds of the different MHD modes become close and the energy can be 
transferred between the different branches of the dispersion relation. The di- 
rection and the effectiveness of the mode transformation depends, among the 
other parameters, on the wave frequency and the attacking angle between the 
wave vector k and the mag netic field (e.g. ICallvll2005l . 120061 : ICallv & Goossensl 
2008). In the two-dimensional case, an approximate for mula fo r the t ransfor- 



mation coefficient from fast to slow mode was derived by ICallvl (|2005T ) for the 
high-frequency waves (above the acoustic cut-off) in a vertical magnetic field: 



exp 



-kir sin 2 ip/ 



dz v\ 



(1) 



where if) is the angle between k and the magnetic field B. According to this 
equation, the fast-to-slow mode transformation is complete (T = 1) for the 

waves with k directed along B. For larger angles ifi the efficiency of the fast- 
to-slow mode transformation rapidly decreases. The higher the frequency of 
waves (implying generally smaller k), the smaller is the cone of ^'s where the 
fast-to-slow mode transformation is effective. 

Numerical simulations of the p and / mode interaction with a vertical mag- 
netic field concentration have demonstrated that, indeed, a significant power of 
the incident modes is converted into the slow high-/3 magneto-acoustic modes 
propagating downwards along the m agnetic field lines and leading to a visi- 
ble power reduction at the surface (|Callv et al.l 11994 ; ICallv &: Bogdanl 119971 ; 
Rosenthal &: Julien 2X100). In addition, these simulations have shown that the 



fraction of the /-mode power transformed to slow modes is sufficient to explain 
its power reduction is observations. On the contrary, the p modes were not 
transformed sufficiently in the vertical magnetic field in order to be explained 
by this mechanism. 
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Theoretical models of the fast-to-slow mode transformation in the inclined 
magnetic field have demonstrated that it can be particularly strong for a nar- 
row range of the magnetic held inc l inations around 20-30 d egrees to the vertical 
(jCrouch fc Canvil2003l : [CaIivll2006l:[Sc"hunker fc Callyll2006h . Similar behavior is 
confirmed also in a 3D analysis bv lCallv fc Goossensl ( 20081 ). The p-mode trans- 
formation at hig her frequencies is foun d to be significantly enhanced by moder- 
ate inclinations ( Crouch Sz Callvll2003l ). The explanation of this effect is offered 
by the ray theory. The fast-mode high-/? waves (analog of p modes) launched 
from their sub-photospheric lower turn ing points reach the tra nsformation layer 
with an angle close to 20-30 degrees (jSchunker Callvi 12006). Thus, for the 
moderately inclined magnetic field th e attack angle ijj in Eq . [J is small and the 
transf o rmation is efficient. Note t hat , ICrouch Sz Callvl ( 20031 ) ; ISchunker &: Calfv 
(12006T ): ICallv & Goossend ((2008) discuss the case where the slow modes pro- 
duced after the transformation pro pagate upwa r ds, w h ile the downward propa - 
gation is observed in simulations bv lCallv et al.l ( 19941 ) ; [Callv Sz Bogdanl (|l997l ); 
Rosenthal fc Julienl (|2000l ). 



The critical role of the magnetic field inclination for the mode transforma- 
tion is confirmed by the numerical simulatio ns of MHD wave propagation in the 
chromosphere bv lCarlsson Sz Bogdanl ( 20061 ). where the authors considered oscil- 
lations with wavelengths comparable to magnetic field scales. At angles smaller 
than 30 degrees, much of the high-/? acoustic-like fast mode power is transformed 
into the low-/? slow modes propagating along the magnetic field lines. At larger 
inclination angles the high-/? fast modes are refracted and reflected and return 
back to the photosphere. Note that the larger inclinations correspond to the re- 
gions where the field strength is smaller and the transformation layer is located 
higher in the atmosphere. When observed at particular height, the interference 
pattern produced by the upward and downward propagating waves in these re- 
gions creates a ring of enh anced power around the magnetic field concentrations 
( Carlsson fe Bogdanl 1200*61 ) which can be an explanation of the acoustic halos 
observed at edges of active regions (Sect. 4). 



6. Local helioseismology in active regions 



Time-distance helioseismology makes use of wave travel times measured for 
wave packet s traveling between various points on the solar surface through 
the interior dDuvall et al l 1 19931: iKosovichevi fl999l. l2002t iKosovichev et al.ll2000l : 
Zhao k Kosovichevl l2003 N ). The variations of these travel times compared to the 



quiet unperturbed atmosphere are assumed to be mainly due to mass flows and 
wave speeds below the surface. The interpretation of the inversion results of 
time-distance helioseismology encountered major critics when applied to mag- 
netic active regions of the Sun. Since magnetic field modifies the wave propa- 
gation speeds in a similar way as the temperature perturbations do, the mag- 
netic an d temperature effects o n the wave travel times are difficult to separate 
(see e.g. iMoradi Sz Callv 20081 ). From observations, it has been demonstrated 
repeatedly that waves travel 20-40 seconds faster deep b elow the sunspots , 



when comp ared to quiet photosphere measurements (e.g. iDuvall et al.l Il993l ; 
lBraun|[l997l ). The explanation of this effect is far from clear. To understand 
the influence of the magnetic field on travel time measurements, the forward 
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modeling of waves in magnetic regions has b e come the prefer r ed approach in 
recent years (iParcheyskv Sz Kosovichevl 120081; lHanasogel 12008 ]_ Cameron et al 



2008UMoradi fc Callv 2008; Moradi et al.ll2008l : iKhomenko et al^ OOS) 



iParchevskv fc Kosovichevl (|2008l ) obtained that only 25% of the observed 
travel time difference can be explained by the magnetic field effects in their 
simulations of waves excited by an isolated source in an inclined magnetic field. 
Other authors report values of the travel time difference s from simulations that 
are similar to the observed ones (ICameron et al.l 120081 : iMoradi fc CalTvl [2008; 
Moradi et aD 120081 : IKhomenko et al.l l2008fi ~ The dependence of the magm- 



tude of the travel tim e diff erence on the magnetic field st rength was shown by 
Cameron etaD (|2008l l and IKhomenko et al.1 (|2008T ). while IMoradi et al.1 (|2008h 



argue travel time differences are largely insensitive to the sunspot structure. 
The frequ ency dependence of t h e travel time measurements reported i n the 
literature (Couvidat et all 120061 ; ICouvidat Raiagurul 120071 : iRajaguru 2008; 
Braun &: Birch 20081 ) is "interpreted as one of the demonstrations of the surface 



magnetic fi eld effects. High frequency waves have, in general, larger travel-time 
differences ( Raiagurul 120081 ) . Qualitatively, such dependence is r eproduced by 



numerical simulations as wel l as the ray theory eikonal approach (jMoradi et al 



20081 : Khomenko et al. 20081 ). The overall agreement between the simulations 



and observations indicates that the observed time-distance helioseismology sig- 
nals in sunspot regions correspond to fast high-/3 MHD waves (modified p and 
/ modes). This conclusion apparently contradicts the chromospheric measure- 
ments, where the correlations between the photospheric and chromospheric sig- 
nals was found and the longitudinal propa gation of the low-/3 s l ow waves allowed 



the be st agreement with the observations (ICenteno et al.ll2 006: Bloom field et al. 
2007al ). In the same vein, ISchunker et"aD (|2005MSchunker Callvi (120061 1 argue 



that the observed dependence of the Doppler signal on the line-of-sight magnetic 
field inclination angle requires an alignment between the magnetic field direc- 
tion and the wave propagation, only possible for low-/3 slow waves. More work 
in needed to clarify these issues including the calculation of the line formation 
heights and sunspot Wilson depression effect. 



7. Conclusions 

Theoretical interpretation of waves observed in sunspots has made big advances 
during the last decades. The influence of the (5 = 1 level and the magnetic 
field inclination on the observed wave properties is being clarified, allowing the 
interpretation of the observed oscillations in term of MHD wave modes. The 
observed wave pattern in the upper photosphere and the chromosphere is found 
to be compatible with the low-/3 acoustic-like slow mode propagation along the 
inclined magnetic field lines. Several physical effects are proposed to be respon- 
sible for the frequency change with height in the umbra from the 3 mHz in the 
photosphere to 5-6 mHz in the chromosphere. The mode transformation at the 
j3 = 1 level is proposed to be responsible (at least in part) for the wave power 
distribution over the active regions, including the wave suppression in the umbra 
and, possibly, enhanced power of acoustic halos. Theoretical interpretation of 
the local helioseismology measurements including magnetic field has been ini- 
tiated. The first simulations of helioseismic waves in sunspots show a general 
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agreement of the wave travel times of the high-/? fast modes with observations. 
Still, the final issue in this field is far from clear and the major developments are 
expected in the future. The 3-dimensional wave propagation and mode trans- 
formation in sunspots are important to analyze from simulations, including the 
chromospheric layers. 
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